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SYITQPSIS 


PmfiMETRIC Brsa?AB3I.IIT OP A PERIODICALLY 
SUPPORTED PIPE CONVEYING PLUID 

This thesis is a study of the parametric instabilities 
of a periodically supported pipe conveying a pulsatile flow. 
The velocity of the fluid is considered to have a harmonic 
fluctuation over and above the mean constant value. This 
work demonstrates that the method of propagation constants 
can be used efficiently for analysis of periodic structures 
vulnerable to flow induced parametric vibrations. In the 
same range of frequencies the number of instability zones 
increases with larger number of spans. The power of the 
method lies in rendering the ccraputational effort independent 
of the number of spans. In addition, zones of instabilities 
associated with all the modes are determined simultaneously. 
Bolotin's method has been used to investigate the bomds for 
the principal primary regions of instability. Sin^e span 
pipe receptances used in the analysis have been obtained 
without taking recourse to mode shape approximation. The 
coupling effect due to fluid inertia has been taken into 
account systematioslly. It is shown that the instability 
zcaies can- be controlled effectively by the use of damped 
dynamic absorbers. 



CHAPTER - I 


lETRODUCTIOU 


1*1 Introduction 

Pipelines conveying fluids at iii^ velocities are 
encountered in various fields of engineering, To name a 
few, such, pipelines exist in the coolant channels of nuclear 
reactors, high pressure boilers, process industries and fuel 
lines of aircrafts and missiles* Due to fluctuations in the 
fluid speed (generated for example by the fluid driving 
pump) such systems are vulnerable to parametric instabilities. 
In severe cases, excessive vibration of the pipelines may 
even cause fatigue cracks and eventual failure of the pipes, 

failure of the coolant channels in a nuclear reactor 
has hazardous consequences. In addition to hazard, the 
hi^ cost of equipment in special applications denands a 
hi^ degree of reliability in the subsystems. Due to these 
reasons, a thorough understanding of the stability domains 
of these pipelines is necessary to ensure proper design. 

With increasin^y hi^ flow speeds required in 
engineering applications, considerable amoiont of research 
work has been done in recent years, in the field of flow 
induced vibrations. A variety of instabilities have been 
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uncovered at hi^ flow speeds that have gone in making this 
subject a theoretician’s delict and a designer’s ni^tmare. 

A fluid flowing (under suitable conditions) in a flexible 
pipe is kno-wn to induce static buckling (divergence instabi- 
lity), combinational resonances, parametric instabilities, 
etc., to name a few. 

Host of the volume of literature is devoted to the 
analysis of a single span pipe. In actual practice however, 
more often we encounter multispanned pipes where analysis 
becomes increasingly difficult with increasing number of spans. 

The present work is a study of the parametric 
instabilities of a unifonn pipe placed on equispaced supports 
and carrying a pulsatile flow. Multispan pipes display more 
number of instability zones in comparison to a sin^e span 
pipe in the same range of frequencies. The method of approach 
used is one that renders the computational effort independent 
of the number of spans. In addition, it enables one to obtain 
the zones of instability, associated with the higher modes of 
a multispan pipe, without involving any further steps in the 
analysis. The objective and scope of this thesis is outlined 
in detail in a later section, 

1. 2 Review of Previous ¥ork 

Since the present work deals with the dynamics of 
pipes conveying fluids and analysis of periodic structures. 
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a brief review of the available literature in each of these 
fields is given in the following sections, 

1.2.1 Dynamics of Pipes Conveying Fluids 

Interest in the study of dynamics of elastic pipes 
conveying fluids was reactivated in 1950 by Ashley and 
Haviland [1], in connection with the vibration problems of 
Trans Arabian pipelines. It was observed that the effect 

of the flow is to damp the transverse vibrations and that the 
flow induced damping (a non viscous effect) decreases the 
response at hi^er flow rates, long [ 2 ] pointed the errors 
in their analysis and showed that damping by internal flow 
occurs only in a certain range of flow velocities. Housner [3] 
reported that at sufficiently high flow velocities the pipe 
may buckle, essentially like a column subjected to axial 
loading. Gregory and Paidoussis [ 4 ] showed for high flow 
rates oscillatory instabilities are possible and that the 
dynamical problan is independent of fluid friction. An 
extensive literature survey in the subsequent developments 
is presented in reference [ 5 ]- 

Chen [ 6 ] studied the problem with the velocity of 
the flowing fluid having a harmonic fluctuation over and above 
the mean constant value. He fo\md that both parametric and 
combinational resonances are possible. He determined the 
instability regions for a simply supported pipe. The 
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governing equation was revised by Paidousals and Issid [7] 
and the regions of parametric instability for single span 
pipe with various end conditions were obtained by the Eitz- 
Galerkin method. 

1.2*2 Analysis of Periodic Structures 

Wave motion in periodic systems in the field of 
crystals, electrical circuits, etc. has been studied for 
nearly 300 years, as pointed out by Brillouin [8] in his 
classical work. The extension of this method for the analysis 
of wave motion in engineering structures (consisting of 
beams, plates, etc.) is a recent phenomenon. 

Mead [9] studied the free wave propagation in 
infinite beams on periodic supports using the propagation 
constants approach. Based on the wave approach Sengupta [10] 
developed a graphical method, to determine the natural 
frequencies of a finite periodic beam. 

Singh aud Mallik [ll] studied the wave propagation 
and vibration response of a periodically supported pipe 
conveying fluid. They showed that due to the presence of 
the Coriolis term in the equation of motion, classical normal 
modes (standing modes) do not exist. Conditions for both 
divergoice instability (buckling) and oscillatory instability 
(resonance) were obtained by the wave approach. 
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A succeeding paper by the same authors [12] was 
devoted to obtaining the instability regions associated 
with parametric excitation of the same system. In this 
work it was shown that all the instability zones are 
insensitive to the fluid mass ratio parameter so that the 
graphical approach is valid for determination of the instabi- 
lity regions. The analysis was based on the governing 
equation as obtained by Chen [6]. Based on the .assumptions 
in [12], the effect of a dynamic absorber on the instability 
zones of a single span pipe was reported by Sin^ and 
Mallik [13] • Mallik and Narayanan [14] extended the analysis 
for a periodically supported pipe with a vibration absorber 
in each span. It was shown that the instability regions can 
be controlled effectively by damped absorbers, 

1. 3 Objective and Scope of the Present Work 

In all these works on parametric excitation of 
periodically supported pipes, it is observed that in the 
initial steps of the wave analysis, certain assumptions are 
incorporated that implicitly ignore the coupling terms. The 
results in these analyses are consequently insensitive to 
the change of the degree of coupling (that is given by the 
fluid mass ratio parameter. ) . Thou^ these assumptions. 


simplify the analysis greatly, this thesis bears it out that 
certain important conclusions reached in these works are 
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erroneous. It is also ranarked that the governing equation 
used in these works is based on that of references [16] and 
[ 15 ], which have been shown to be incorrect [7]. 

The present thesis is an attempt to rectify these 
shortcomings in the previous analyses by a rigorous 
application of the propagation constants approach. Moreover 
the correct equation of motion which has non-constant 
coefficients is used. The coxipling effects have been 
retained at all stages of the analysis, and it is shown 
under what restricted situations and stages of analysis are 
the conclusions of the above mentioned references justif ia,ble. 

As mentioned before, the ccmputational effort, in the' 
methodology used, is independent of the number of spans. 
Computations have been mdertaken only for the case of one 
span and two span pipes for the sake of brevity and for 
comparison with the results of previous works. 

As long as the periodicity in the structure is 
maintained, the method can be used for any tsrpe of support 
conditions. For the sake of simplicity it is assimied in 
this work that the pipe is placed on line supports (i.e. no 
axial force and no torsional restraint is produced by the 
supports). Due to complexity of the governing equation, no 
mode shape approximations have been resorted to. The 
transverse displacement has been expressed as a truncated 
power series of suitable order so as to obtain the desired 


accuracy 
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The effect of attaching a dynamic absorber at 
the midspan of each bay, is also investigated. It is 

demonstrated that, by suitably changing the absorber 
parameters, the instability zones (in excitation parameter/ 
frequency space) can be effectively controlled. The matter 
has not been explored beyond this stage since optimisation 
of the absorber parameters would be a separate topic in 
itself. The analysis has been restricted only to the 
determination of the principal primary instability regions. 

The present thesis may be broadly divided into 
three sections : 

i) Determination of propagation constants for an 
infinite periodic pipe. 

ii) Determination of the zones of parametric 

instabilities of a finite periodically supported 
pipe. 

iii) Investigation of the effect of dynamic absorbers 
on the instability zones of these pipes. 


« 



i 




CHiPTEil - II 


PAKiMEIHIC KS TABU. in OF A PERIOD IC ALB Y SUPPORTED P IP E 

COWEYIHG FLUID 

2* 1 Introduction 

In this chapter we shall discuss the means of 
obtaining the instalility bounds of a parametrically excited 
pipe placed on equidistant line supports. For a pipe of N 
spans we have adopted the convention of naming the supports 
as 0, 1, 2,..., R as shown in Fig. (2.1). In section 2.2 the 
governing equation for such a pipe carrying time dependent 
flow is discussed. We assume the flow to have a harmonic 
fluctuation over and above a mean velocity. In section 2,3 
the goveming equation is reduced to a form amenable for 
further analysis, especially for the determination of regions 
of instabilities. A periodically supported pipe of infinite 
number of spans can be considered to be composed of identical 
members, each of one span length, connected together to the 
adjacent members in an identical manner. In section 2.4 
using this spatial periodicity, a functional relationship of 
the propagation constants in terms of the single span recep- 
tances is obtained, so as to satisfy the intermediate support 
constraints. Section 2.5 is devoted to develop the equation 
for obtaining the frequencies bounding the instability zones, 
by application of the end support constraints for an H-span 
pipe. In section 2.6 an outline of the method to obtain the 
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receptances for a single span pipe is presented. A particula 
case of parameter values for which the analysis is simplified 
is discussed in section 2.7. lastly in section 2.8 the 
discussion of results and computations is undertaken. 

2.2 Equations of Motion 

The equation for transverse motion of a pipe carrying 
a time dependent flow, in the absence of external loading is [ 


I li-w— + + (m„+m )g’ ^ + 

3Z 31 ; d'5. ^ 3X 

[m^-U -■ T + pA (1 - 2 vf) - ((m^+ni )g’ ~ m^ ]”ir 

i 9 2C 

+ 2 u i|i- + 0 ^ + (m.+m ) -^ = 0, 

^ 9x3t ° ^ 9T 

0 <_ X < L , ...(2.01) 


where, 

T is the independent time variable, 

X is the independent length coordinate along the 
pipe axis, 

1 is the unsupported pipe length between two adjacent 
supports , 

w is the transverse pipe displacement and is a function 
of both X and %, 

I is they(^cross sectional moment of inertia. 
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is the viscoelastic coefficieat of internal dissipation, 

E is the elastic modnlns, 

is the mass of the fluid per unit length, 

is the mass of the pipe per unit laigth, 

g’ is the axial component of the gravitational acceleration. 

U is the time dependent bulk velocity of the fluid, 

I is the externally applied tension, 

p is the fluid pressure, 

A is the internal pipe cross sectional area, 

V is Poissons ratio, 

f=0 with no axial constraint at the supports, 
f=l with axial constraint at the supports, and 
c is the coefficient of viscous damping applicable to 
transverse motion of the pipe. 

Equation (2.01) was derived with the following 
assumptions % 

(1) The amplitude of transverse oscillations is small 
with respect to the length of the vibrating pipe. 

(2) The fluid is inviscid and incompressible and has a 
plug flow profile. 

(3) The flow is iminfluenced by the transverse motion 
of the pipe. 
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To simplify the analysis we retain only certain 
terms and make the following additional ass’umptions z 

(1) low damping in pipe material, implying B* = 0. 

(2) Pipe is li or i a on tally placed, for which case g’ = 0, 

(3) No azial tension, or T = 0, 

(4) Negligible effect of pressurisation, i.e. p = 0, 

(5) External viscous damping is very small, or c = 0. 


With those assumptions, the governing equation (2.01) 
reduces to. 


El + m^ w . r.. -N 3U-^ a^ w . _ ^ 

g-4 f + Lm.(L-X} — J — ^ + 2m^ U 

^ 3 X^ ^ ^ ax 3 i; 




w 


+ (nif + m ) 2 


0, 0 £ i 1 1 . . . (2.02) 


The non-dimensional form of equation (2.02) is. 


w . „2 3^w 


^ Ud-x) ^ + 2y„ ^ H. = 0, 


3X 


ax 


3X 


a w , a‘~w 
at' 


0 


X < 1 


...(2.03) 


where the following non-dimensional parameters are used: 


(1) Transverse displacement w = w/L, 


1 ^ m , 


(2) Eluid velocity u = U 
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( 5 ) 

(4) 

(5) 


Axial coordinate x = x/L, 


Time t = t [ 


M ii/2 


(m^ + Dip)!' 


Fluid mass ratio parameter y = [ 


m. 


f p 


, 1/2 


In equation (2«03) the leading term represents the flexural 
force and the fifth term represents the total inertia force 
due to transverse motion of the pipe and fluid. The second 
term arises due to reaction on the pipe hy the fluid as it 
moves in a curved path due to pipe flexure (centrifugal 
effect). The third term arises due to the time dependent 
nature of the flow. Also as the pipe is undergoing transverse 
vibration^ at each point the slope of the pipe changes with 
time (and hence -the fluid is forced to rotate at each point) 
causing a coriolis force to act on the pipe represented by 
the fourth term. This coriolis term causes coupling between 
the Spatial modes of the pipe, the implications of which, will 
be evident later. It may be noted that the second term is 
equivalent to an axial compressive load on the pipe and for a 
suitably high magnitude of u causes static buckling [llj. 

In this work, a pulsatile flow has been modelled as a 
non- zero quiscent flow carrying a harmonic rider which is 
represented thus, 

= Uq (1 + 6 cos nt) 


u 


. ..( 2 . 04 ) 
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where is the quisceat fluid velocity and 6 is the 

excitation parameter, and is the non-dimensional 
frequency given hy with u rad/time as the dimen- 

sional frequency. 

Substituting equation (2.04) in equation (2.03) 
one gets, 

4 _ 2 

+ [^ 0^1 + 6 cos fit) - j Uq 6 fi(l-x)sin fit j 

3 ^ ^ 3X 

2 2 

+ 2yu„(1+6gos fi t) '-j ^ 4 : + = 0, 0 < X < L 

o gx 3 “C 

...(2.05) 

The coefficients of various terms in equation (2.05) are 
time varying and hence the system is said to be parametrically 
excited. With constant parameters (i.e. coefficients) the 
system would have had instabilities at certain discrete 
frequencies. The effect of parametric excitation (due to 
time dependent velocity in our case) manifests itself by 
causing the system to have unbounded solutions (instabilities) 
at certain discretely spaced frequency bands for a given 
value of the excitation parameter, 6 . The unstable regions 
are normally plotted in the 5^-6 space. 



15 


2,3 Determination of Regions of Parametric Instabilities 


Ihe regions of primary instability of equation ( 2 <. 05 ) 

Pt 

are bounded by periodic solutions of period 2IE where T = 
Hence, on the bomdary, the transverse displacement can be 
expressed as [16], 


w(x,t) = 


I (Xg(x)sin(-|^) + yg.(x)cos(-|^)) 

... (2,06) 


X=1,3,5,.. 


The region of principal primary instability can be obtained 
quite accurately by truncating the series in equation (2.06) 
at K = 1 itself [l6]. Hence one can write 


w(x,t) = X(x) sin(~^) + Y(x) cos( 


r Qt v 




...(2.07) 


where, for convenience we dispense off the subscripts on 
Z(x) and Y(x). In this approximation we are neglecting the 
instability zones arising due to K = 3»5>*..« This step is 
justifiable since the higher instability zones are very narrow 
and are therefore of- negligible interest. Also, by retaining 
the terms with K > 1 in the equation (2.06) for w(x,t), 
the principal primary Instability zones are progressively 
modified, but the change whatsoever is small and may be 
ne^ected. Substituting equation (2,07) in equation (2.05) 
and equating coefficients of Sin (“^) and Cos(— separately 
to zero, we get. 
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+ u (1-6+ Y -) 7-2 

dx^ dx 


7 

4 


X" ^ 4 iL [(i_x) ifr 


dx‘ 


_ (1 . 2x din 

6^ dx-> 


= 0, 0 £ X £ 1, 


. . . ( 2. 08a) 


d^ 

dx'^ 


^ (1+6+ ) 


yu^bQ 


dh 

dx^ 
d^X 


i n^Y 

4 


n + 2x ^n 
+ 6^ & 


0, 0 < X £ 1, 


,..(2.086) 


These are two fourth order coupled ordinary differential 
equations with non-constant coefficients. The coupling 
between the equations of X(x) and Y(x) originates from -the 
terms arising due to unsteady flow velocity and the coriolis 
acceleration. For y - 0, equations (2.08a and b) are decoupled 
and analysis becomes simpler. This case is discussed in 
section 2.7. 

It is shown in reference [5] that the influence of 
neglecting the coupling terms (i.e. making the y = 0 approxi- 
mation for the situations in which y is actually non-zero) 
does not change the instability zones significantly. However 
due to seme inconsistencies in the analysis presented there 
and in view of the improved governing equation (2.01) which 
is different from that used in reference [5], the effect of y 
on the stability zones needs fresh investigation. With this 
aim the coupling terms are retained for further analysis. 



17 


For a given set of parameter values the instahility 
hounding frequencies are those at ishich equation (2.08a and h) 
are satisfied alcngwith the boijndary conditions of the pipe. 
For example in the case of a single span pipe placed on line 
supports at both ends the relevant boundary conditions at 
X = 0 and x = 1 are, ' 


w(o,t) = 

X(o) 

sin(— 1^) + Y(o) cos(~^) = 0 

(2.09a) 

w(l,t) = 

X(l) 

s±n(|i) + Y(l) cosCf^) = 0 

(2.09b) 

^(0,t) = 

X (o) 

sinC”^) + Y (o) cos(^) = 0 

(2.09c) 


2 

■^(l,t) = x"(l) sin(|^) +Y“ (1) cos(-~) = 0 (2.09d) 

3x 

In terms of X and Y equations (2. 9a,b,c,d) may be written as 


X(o) = 

Y(o) = 

0 

(2.09e) 

X"(o) = 

Y*'(o) = 

0 

(2.09f) 

X(l) = 

Y(l) = 

0 

(2.09g) 

X”(l) = 

Y>'(1) = 

0 

(2.09h) 


Thus for a single span, the systeu of equations (2.08a,b) 
consists of 8 unknowns and 8 boundary conditions ( 2.09e,f , g,h) . 
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By solving ior the eigen values of the coupled equations 
(2.08a,b) with the appropriate boundary conditions, we obtain 
the frequencies at which a non-trivial solution exists. These 
frequencies then corresponds to the bomding frequencies of 
the instability zones for a single span. Extension of this 
method (referred to as the direct method) for a pipe consist- 
ing of a large number of spans is mweildy as 8 boundary 
conditions are to be satisfied at every intermediate support 
over and above satisfying four boundary conditions at each 
end support. Thus for an h-span pipe, the frequency deter- 
minant (zeroes of which gives the eigen values) is of the order 
8IJ X aj. The direct approach becomes unwieldy due to the 
appearance of numerous intermediate support conditions, for 
large number of spans. 

However, if an N-span pipe is periodic in nature 
i.e. consists of identical elements connected together in 
an identical fashion, then the intermediate conditions between 
any two elanents are identical in nature. In such a case, all 
the information contained in the many intermediate supports, 
may be substituted by obtaining a set of propagation constants 
(to be explained la.ter) by satisfying all the constraints at 
any one of the intermediate supports. By this approach the 
computational effort is rendered independent of the number 
of spans. This method using the concept of propagation 
constants is outlined in the subsequent section. 



Direct ion dt flow 



I 


F»g,2-2fb) Two adjacent periodic elements the infinite pipe 
at the I N ) th support 







2.4 Determination of Propagation Constants for an 
Infinite Periodically Supported Pipe 

Let us consider a pipe, infinite in extent, placed 
on equispaced line supports as shown in Pig. (2.2a) and 
carrying a flow described by equation (2.04). Equations (2.08a, 
b) then apply to each span of this pipe. We shall satisfy the 
slope continuity at say the Nth junction in the following 
manner. Consider the two spans AB and CD, adjacent to the Nth 
support as shown in Pig. (2.2b), Let be the non- 

dim ensionalised moments (Mj^ = , where is the dimen- 

sional moment) acting on the two ends of the span AB. ¥e have 
satisfied moment continuity relation at the Nth support by 
assuming to act on both B and C. let 0g and represent the 
slope of the pipe at endsB and C respectively. Por the system 
response i w(x,t) given by equation (2.07) we may consider the 
slope and movement at anjriocation to be similar harmonic time 
functions of frequency ( *^ ) . Hence we designate the following 

coefficients (unknowns); 0^ i q= 2 :,y 

111 11 

such that ; 

( — )b = sin ( y ) + 


. . . ( 2.0 9a) 
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(m) - 

®0 = ®0x + % OOBif) 

...( 2 . 09 b) 


= sin(fS) + oos(|i) 

. . . (2.10a) 


= oos(|i) 

... (2,10b) 


= nf-l) eiuif) + cosif) 

...(2.10c) 


Using the influence coefficients or receptance method, the 
slope at B may he expressed as a linear combination of the 
end moments acting on the span AB. Due to coupling in the 
governing equation (2.0ea,b) a sinusoidally varying end 
moment in addition to exciting a sinusoidal response also 
generates a cosine response. To this end we may distinguish 
the receptances for a single span pipe in two categories. 

(1) Direct receptances : 

Response at A and B to a imit moment acting at one 
of the ends which is in phase with the exciting moment, 

(2) Cross receptances ; 

The response at A and B which is out of phase 
(by + with the exciting moment acting at either of 


the two ends 
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With this in mind may he expressed, in general, 

as a linear combination of the four tjmes of moments 

^x^ ' ^^y^ My^^^ acting on pipe element AB, as 

below s 


9. 


Bx 


pS ^ - pS Mf ^ - pS - pS 


, . . (2,11a) 


©-^ = 
By 


Pll ^^X ^ Pll \ + PlO ^ PlO 


'11 


...(2. lib) 


where for example P^q is defined as the amplitude of the 
sinusoidally varying slope response at the right end (1) of 
the span generated by a unit cosinusoidal moment acting on 
the left end (o) of the span. The remaining receptances are 
defined and scripted similarly. The direct receptances are 
identified as and p?^ , i = 0,1; j = 0,1 and the cross 

X J X J 

receptances correspond to the type P?^^ and p^^., i = 0,1 
and j = 0,1. The receptances p's are properties of a sin^e 
span pipe as governed by equation (2,08a, b) and are functions 
of the non-dimensional frequency fl and other systm para- 
meters only. The method for ccmputing the receptances is 
outlined in section 2,7. Since the spans AB and CD have 
identical physical parameters a similar set of receptances 
may be used to express 9^^ and ©^y in terms of the end 
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moments acting on span CD, Hence 


® 0 x = Poo 4 “’-^ Poo 4 ^’^+ PS 4 “'*'^^+ pS 4 ”"^^' ...( 2 . 12 a) 


O- 


!y = Poo Poo, Pm ...(2.12b) 


’ 00 . 


'01 “x 


01 “y 


It is remarked that 4- due to the i 

directional assymmetry imposed hy the flow and the consequent 
appearance of terms depending on u, ■& in the governing 
quation (2.03). However a restricted sort of equivalence 
^or hold het^reen certain cross receptances and this is 
’....icussed later section 2.7. 

Dor the condition of slope continuity = ©q» to he 
true for all time, the sine and cosine components of the two 
hmetions must necessarily he equal, i.e. , 


Sx “ 

®Cx 

... ( 2 . 13 a) 

®By = 

Cy 

...( 2 . 13 b) 


Substituting the expressions for ©^^ and ©^^ from equations 
(2,11a) and (2.12h) in (2.13a) and subsequently rearranging 
the terms one obtains 

XX j^(H-l) /.xx.xxn j.(N)_ gxx j.(N+l) „ gXy (H+l) 

•^10 +^Pll POO^ POl ^^x ■" POl 


• • • ( i2* 14*0'} 


+ (pS " pS) 
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and similarly from eq^uations (2.11b), (2.12b) and (2.15b) 
WQ get, 




11 


00 ^ 


’01 ""x 


+ (pgg - pg) pSm(®-1> 



M 


(N+1) 

y 


...(2.14b) 


Due to spatial periodicity [9], a relationship which holds 
between two quantiuies (e.g. slope or moment) at one pair 
of adjacent supports must necessarily also hold true at all 
other pairs of adjacent supports as well. Thus equations 
(2.14a) and (2.14b) can be written for any three consecutive 
supports. Solution of these coupled difference equation 
may be written in the following form 



h(2) 

,(H) 

X 



^(N-l) 

^ X 

X 



yi(M) 

_y 

m(H+1) 

. y .. . 



•••• >-i) 





... (2.15a) 


J-P' 


...(2.15b) 


where ii, is termed as the propagaticai constant [9]. 

In general, p is a complex number where the real part 
signifies the logarithmic decrement per span and the imaginary 
component signifies the phase change per span. 
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Substituting the moments at the (H+l)th and (N“l)th supports 
in terms of moments at the (N)th support in equation (2.14a,b) 
we obtain 


(pS-PS)-(e^) - 

. . . ( 2. 16a) 



2.16b) 

Equations (2.16a,b) are two homogeneous linear equations in 

and renresenting the condition of slope continuity. 

X y 

For a non- trivial solution the determinant of the 
coefficients of equations (2.16a,b) must vanish. Hence we 
obtain the condition, 



.(2.17 a) 
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and 


jyj(N) 

WJ 


M 


r(^i) = 


X 




00 ^ 11 ' 


(P5(e'^) ^(Pgg-Pg)(e>^)-Pg) 
(pgf(e‘^)^(pgg-pg)(e>*)-pfg) 

...(2.17b) 


where, r is defined as the ratio of the out of phase moments 
at each support. 

Equation (2.17a) is a fourth order polynomial in ( e^) having 

i^i 

real coefficients. The four roots of this equation e , i=l,4 
correspond to the four values of the propagation constant. 
Moment continuity has already been incorporated as discussed 
earlier. Also the constraint of zero displacement at the 
supports has been taken care of by using appropriate receptan- 
ces (section 2.7). Since there exist four roots e , i=l,4,for 
each p. > we obtain the cigen-'Vectors :• and sn r^ 


such that 


(N) . (H) 


/n' 

y,i^ x,i i " i' 


= r. = r(Pj^}, i = 1,4 ...(2,17c) 

The moment components at the (]S[)'th support may be expressed in 
general as a line.:r combination of the eigen-vectors. Hence 

. . . ( 2.18a) 


(N) _ 

l^(b) 

+ M 

(N) 

+ + 

j^Ch) 

x “ 

X,1 

x,2 

x,3 

x,4 

(H) _ 


+ M 

(H) 

+ + 


y 

y,i 

y»2 

y»5 

y»4 

= 

ri-M 

(H) 

X,1 

+ r 

2 x,2 

3 


4 x,4 

..,( 2.18b) 
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/ |,T N 

where, a,refoi:ir unkiiown momeats of the same phase 

at the Nth support. 

Using equation (2,15a and h) we may relate the moments at 
the (N+l)th and (N)th support as, 


M 

M 


(N+1) 

x.i 

TST” 

x,i 




r 1 , 4 ^ N"0, ^1 , ^ 2 , • , • 


. ..(2.18d) 


Therefore, in general at the (N+l)th support 


M, 


X 




...{2,19a) 


and 

jy[(N+l) 










...(2.igb) 


Thus we have satisfied the intermediate constraints by 
prescribing 

(i) the propagation constants which relate the mcmaits at 
different supports (equation (2,17a)), 

(ii) the inter-relation between the out of phase moments at 
a support (equation (2,18c)), 



Direct ioR of flow 
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For complete solution there remain only four unknowns 
i. e. ^x,2’ 3’ 4 determined. 

As shown in Fig, (2,3) if the infinite pipe under considera- 
tion were constrained to oscillate in such a fashion that the 
moments at the (0)th, (+11) th, (+aj)th, ... supports were zero, 
the system would be equivalent to a periodic extetision of a 
finite H-span pipe having zero end mcments. Hence the method 
of propagation constants developed in this section may use- 
fully be applied to obtain the bounding frequencies for a 
finite N-span pipe by imposing the condition of zero end 


moments. 
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2.5 Determination of Primary Instability Boundaries 


In this section, we shall discuss the method of 
obtaining the principal primary instability boundaries of 
an D-span periodically supported pipe. In other words, we 
will seek those frequencies JCl- , 3 = 1,2,3.*.> (termed 

J 

bounding frequencies) which will generate such propagation 
constants i = 1,4 that cause the moments at the end 
supports to vanish. 

Let denote the total moments at the 

extreme supports (see Dig. 2*5)* Then using equations (2.18a, 
b,c) and (2.19a,b) we may write for as, 

sin (^) + cos i^) 

- (M^ + M2 + + M^) sin (-l^) 

+ (M^-*r3_ + M2*r2 + ^3*^3 + ^4*^4) ^2^^ 

...( 2 . 20 ) 


where for simplicity we denote 



M 


( 0 ) 

x,i • 


At the Nth support 


M 


(N) 

x,i 



p. N 

M^*(e ^) 


(N) ^ M^^^-(e^^)^ = r.-M.*(e ^) 

y,i y,i ^ 


...( 2 . 21 a) 


i=l ,4 ..( 2 . 21 b) 



Hence, 

= H») Sin (|S) H- m'**) cos (|i) 

N li N H 

= [M^-(e ) +M2 *(e^) +M^*(e^) + 

^ o+ 

+ M4 -(e ^ ]sln (f6) 

+ [ r^*M^*(e ) + r2*M2*(e ) + r^*M^*Ce 

+ r^-M^-(s4) j oos (ft) 

...( 2 . 22 ) 

For zero end moments the coefficients of sine and cosine 
terms are separately equal to zero in equations (2*20)and 
(2.22). Therefore, 

Ml + + M 4 = 0 ...(2.25a) 

^1*^1 ^2*^^ 2 ^3*^5 ^4*^4 " ^ ...(2.25b) 

M N p H 

(e ) + (e ^) •M2 + (e ^) ‘M^ + (e ^) = 0 

. . . ( 2. 24a) 

P-j H Pp N p.^ H p* 

ri-(e ) -M^ + r2*(e ^ •M 2 + r^‘(e ^) -M^ + r4-(e ^)-M4 = ( 

...(2.24b) 


For non- trivial solution, the determinant of the above system 
of equations must be equal to zero i, e. , 
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!'(«) = 


(e 


p.. N 


Pp H 

(e 2) 


p, N 

(e 5) 


li, H 

(e 


t-1 - P^ N P^ IT P . N 

r-i'Ce ■^) r2*(e r^*(e 

...(2.25) 


= 0 


The frequencies j = 1,2. ..., where = 0 

3 

correspond to the bounding frequencies of the instability 
zones. 


The methodology used in the analysis may be summarized 
as follows. At the boundaries of the instability zones the 
system is governed by periodic solutions. For principal primary 
instability the solution is approximated as, 

w(x,t) = X(x) sin i-^) + X(x) cos 


For such a form of the solution, we obtain receptances for a 
sin^e span. Propagation constants are obtained by making use 
of the receptances and the constraints at the intermediate 
supports. By requiring the end moments to vanish we arrive at 
a constraint in terms of the propagation constants. Those 
frequencies at which this constraint is satisfied, correspond 
to the frequencies bounding the instability zones. 



33 




2.6 Determination of jR,eceptan.ces 


In the previous analysis, sixteen different 
receptances are usea to relate the slope and the internal 
moments at idle supports. These receptances are obtained by 
solving for the response of a single span pipe subjected to 
four different combinations of harmonic end moments. 

Consider a single span pipe as shown in Fig. (2,4) 
governed by the equations (2.08a and b). The displacement 
boundary conditions may be expressed as (see equation 2,08c), 

X(0) = y(0) = 2 ( 1 ) = Y(l) = 0 ...(2,26) 

Let the syst®a bo acted upon by a unit sinusoidal 
moment at end A, Thus 

= 1 • sin { -Si ) 

° 

The dynamic response of the pipe, which is in linear proportion 
to the magnitude of the exciting moment, will be composed of an 
in phase and a quadrature component. The amplitude of the 
slope response at the two ends A and B, which is varying in 
phase with the sinusoidal moment at end A corresponds to the 
receptances Similarly the receptances and 
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are the quadrature ccmponents of the slope response at ends 
A and B respectively. Hence, one can write 



Sin ( f ) + pgg cos ( f } 


( ~ ) 
3X 


f ^ ^ \ f S^t \ , / 3 Y \ ^ „ / Pt \ 

^ 3 X ^B ^2 ^ ^ 3x ^B ^ 2 ^ 


= Sin ( cos ( ^ ) 


The ei^t Boundary conditions for this case are 


X(0) 

= Y(0) 

II 

H 

II 

H 

11 

O 

... (2. 26) 

X''(0) 

= -“x 

= -1 

. . . ( 2. 27a) 

Y'’(0) 

= - 

*7 

= 0 

... (2. 27b) 

X»(l) 

= - M® 

X 

= 0 

. . . ( 2. 28a) 

Y’'(l) 

= - 

y 

= 0 

... (2. 28b) 


In the subsequent analysis we obtain a solution for this 
systm of equations by the power series method. 
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L et , 

K 

X(x) = -I (1-x)^ ...(2.29a) 

in=0 

K 

Y(x) = -I (l-x)“ ...(2.291 d) 

in=0 

¥e substitute equations (2.29a and b) in (2.08a and b) and 
equate coefficients of like powers of (1-x) to zero, Ibis 
translates the differential equations into the following 
recursive relations. 
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Sub stitu ting equations (2,29a and b) in. equations (2,26, 
2,27a and b, 2,28a and b) we may rewrite ube boundary 
conditions as. 


E K 


1 


=: 

-I 


So = ^0 = 0 

... (2.31) 

m=0 



m=0 




K 







1 

m(m- 

-1) 



= -1 

. . . (2.32a) 

m=2 







K 







1 

m=2 

m(m~l) 

\ ~ 

-«x 

= 0 

... (2.32b) 

2 g2 


- I' 

([B 

X 

0 


... (2.33a) 

2 h2 


- I' 

ii 

0 


, « • ( 2, 33b ) 


The number of unimowis (g^, and h^, h^,....,hg) 

totals to 2 (K+l), The set of two recursive relations 
(^ ” ^ ♦ 1» 2, • « » , (E“4') ) gsn erat e 2 (£*“3) equati ons and with 
8 botindary conditions the total number of equations is 
2 (E+1), and is equal to the number of unknowns. 

Thus the problaa has been reduced to finding a solution 
for the set of linear equations for which standard algorithms 
are available. Once the coefficients h^j m=:0,,.,.X have 

been computed we obtain the receptances from the following 
relations % 
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E K 

^00 ~ ^3x^x=0 ~ ^ ™ = - I m ^ 

■m — T X=0 -, 


• • • ( 2 « 3I3,) 

iL 

oyx _ „ 

Poo - ^x^x =0 = “ Z m li^ ...( 2 , 3113 ) 

m=l 




XX 

10 


yx 

10 


“ §1 



Case II : 

let the end moments he 


...(2.32a) 
... (2.3213) 


= 1 . cos ( ) 

= 0 


...(2.33a) 

...( 2 . 3313 ) 


For such a case we obtain 


Ppg = X'(0) 


„yy 

Poo 


Y«(0) 


jXy 

’10 


Z'(l) 


,yy 

^10 


the following receptances 

...(2.34a) 

...(2.34b) 
. . . ( 2. 34c) 


I’(l) 


..(2.34d) 



Following the method outlined for the earlier cases 


with M =0 and = 1-sin ( ^ ) we obtain 
and with 


^01’ 


Pi 1 f 


rA 


rB 


= 0 and = l*cos ( ) we obtain 


11 

Poi» 


oxy 

>^11* 


nJJ 

^11” 


11 ’ 

xy 

01 ’ 


It may be shotm using reciprocity relations for dynamic 
receptances that. 


pXX 

>"10 

pXX 

- " Poi 

...(2.35a) 

fyx 

f^oo 

r^y 

'•00 

...(2.35b) 

pyx 

^10 

= - pjl 

...(2.35c) 

pyx 

Poi 

= - pS 

...(2.35d) 

„yx 

^"11 

pxy 

^'11 

...(2.35e) 

pyy 

‘ao 

- - r/y 
>^01 

...(2.35f) 


2,7 Determination of Instability Bounds for the 
Uncoupled Case 


It has been mentioned in section 2.3,'fch-a-t if the 
terms with the fluid mass ratio parameter, y are ne^ected 
in equations (2.08a and b) , the method for determining the 
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instability regions is considerably simplified. Sncb an 
approximation is valid when liiefluid mass ratio, y is much 
less with respect to the buckling ratio as defined below 


where 0 


Y << 0 < 1 


faf Pq) 


and Y 




] 


1/2 


In effect, by reducing y = 0, we have neglected the terms 
arising due to 

(i) coriolis acceleration and, 

(ii) unsteady flow velocity with respect to the term 

2 a^w 

in eqn.(2.03)(u ) representing centrifugal effects. 

3x 

It may be noted that the constraint on 0 (the ratio between 
the equivalent axial load due to fluid flow with the first 
Euler buckling load) is less stringent. (All that is 
required is that 0 < 1). 


Hence substituting y = 0 in equations (2.08a and b) we 
obtain, 


^ + u^ (1 - 6 + 

dx^ ° 


— 2 + u (1 + o + 
dx^ ° 



...(2.36a) 


0 


...(2.36b) 
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Equations (2.36a and b) are two uncoupled equations each 
represaiting simple transverse vibration of a beam with 
different axial loads. The frequency i ^ ) oan be considered 
as the natural frequency of these beam columns. The 
calculations of receptances becomes easy since equations 
( 2 . 36 a and b) may be solved analytically. The cross recep- 
tances for this case will be aero and with reference to 
Pig. ( 2 . 4 ), the pipe slope at support B will be (see equations 
( 2. 11a and b) ) ^ 


M 

11 


+ -U 

X 

■ pS 

...(2.37a) 

II 

>> 

CD 

• p'g 

y 

• Pio 

...(2.37b) 

where, 

-L J 

= = 

0 i = 

0,li j = 0,1. 



By applying the condition for slope continuity in conjunction 
with equation (2.36a) we obtain a quadratic equation whose 
two roots e^^ yield the propagation constants for 2;(x). 

In a similai’ manner equating 0^^ = we obtain another pair 

of propagation constants pertaining to y(x) or the 

cosine component. For comparison with y ^ 0 case it may be 

noted that in section 2.4 we had obtained a single fourth 
order polynomial in ( e^) which yielded the four propagation 
constants. The fourth order polyncmial of the propagation 
constants degenerates into two independent second order 
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polynomials in (e*^) for the .y = 0 case. As in the case 
of a harmonic wave in an infinite periodic beam [9] in each 
set (p.^, P2) P-^) » 'the two propagation constants 

are of opposite sign. (P-p = “ Pi P4 = ” P3^ 


To obtain the instability bounds for an h-span 

pipe we apply the zero end moments condition, Thus if 

( 0 ) ( 0 ) 

1 * 2 nnknom moments corresponding to p^, 

P2 respectively, at the (0)th support 


( 0 ) 

m'^ 

y»2 


M 


( 0 ) 


sin ( ■— ) + cos ( -fi) 

■JO) ^^( 0 ) 


K"i ^x^2^ cos(-f ) 


‘y,l • *^y,2^ 

...( 2 . 38 ) 


Por M 


( 0 ) 


0 we obtain 


+ «x?2 


«y?l + «y?2 


= 0 


0 


...( 2 . 39 ) 


...(2.40) 


At the KTth support in a like manner we obtain 

® 


,«> , «»> 


. . . ( 2. 41a) 


0 


...(2.41b) 
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Due to decoupling we may t-rrite the above eauations as 


(e'^2/ = 0 


y.2 


...(2.42) 


...(2.43) 


Thus from equations (2.39) and (2.42) we obtain 


ti. N 

(e b 


tl, U 

(e 2) 


M 


( 0 ) 

x,! 


M, 


r(0) 

x,2 


f ={0} ...(2.44) 


For a non- trivial solution of equation (2.44) 


1^1 N M-o w 

(( e ^ ) - ie ) 


0 


...(2.45) 


Sengupta [3.0] has shown that for determination of natural 
frequencies only the propagating waves (Ee(ia^) = Ile(p. 2 ) = 0) 
need to be considered. Equation (2.45) can be simplified to 

Re (P-^^) = 0 - ...(2.46a) 

Im (p3_) = m = 0,1, 2,..., IT. ...(2.46b) 

In a similar manner for the cosine component we obtain from 
equation (2.40 and 2.43) 
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Im (ti,) = ^ , m = 0,1,2,..., ...(2.48b) 

J 1 '^ 

Frequencies vjithin the propagation bands where 
conditions (2.46a,b and 2.48a,b) are fulfilled are the 
desired primary instability bounds. For each instability 
Zone one bounding frequency is obtained from equation (2.46a 
and b) and the other from equation (2,48a and b). Equation 
(2.47b) (and also equation (2.48b)) can be solved by a simple 
graphical procedure presented in reference [lO]. First a 
curve of Im(|i' 2 ^) versus n is drawn for 'che propagation band. 
Then by dividixig the range of Im(|i^) (0 to it) into N equal 
divisions, one obtains thenatural frequencies in each band. 

As explained in references [9] [l^,the highest frequency 

in each band, as obtained by this method has to be discarded 
and thus one gets N natural frequencies from every propagation 
band. A similar procedure followed with the curve Im(p.^) 

Q yields the second set of U natural frequencies. 


versus 
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In the following exialysis it is shown that for the 
case, Y=^ method using the sero crossings of PCfi) 0 'ut~ 

lined in section 2.4 is equivalent to the graphical approach, 
por Y = 0, equations (2.08a and h) are uncoupled therefore 


pxy 

1^1 Up 

let e , e 




= 0 


i = 0,1 , 


0,1 


jj, [I 

he one pair of conjugate roots, and e e 


he the other pair. ¥e redefine the ratio r^^^ such that the 
cross receptances appear in the numerator and the direct 
receptances in the denominator (for proper limiting). 

Let 


r ' 

i 


M 


M 


(^) 

m 

x,i 


(p 


,yy 

00 


pS) -Pio 5 


for i = 1, 2. 


...(2.49a) 


and 


r I 
1 


M 


M 


(N) 

^,i 

ity 




- p??)-(e% - pS? ) 


'01 




(Poi*(e ) + (Poo •" Pll^ ■^“‘" 10 ^ 


XX > 


U. 


,xx 


for i = 3,4. 


...(2.49h) 


Ihe condition of zero end moments equations (2, 23a,h and 2. 24a,h) 
in terms of r^, i = 1,4 becomes 
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x,x x^d ^ y,:5 4 7*4 


...(2.50a) 


(/ . (e“ '^2).^(0) ^ ,,.^(0)./ ^,.^^0),-/ ^4 . 0 


‘y»4 

..(2,50t)) 


ri4°{ + r- -H^ + m(°) + « 


. ..(2.51a) 




W^i, 


Wli. 


N^x, 


...( 2 . 516 ) 


■wliich may be written in matrix form 


1 

1 

r ’ 

3 

r ' 

? 


4?i 

. 

e 

e 

Up 

r’-e 

Up^ 

r’-e ^ 

4 

*< 

^4^2 

X,2 ., 

> 

r£ 

r ' 

2 

1 

1 


y,3 

Nti 

ri*e ■*■ 

r^' e 

Up, 

e 

UP^ 

e 


y, 4 : 

— -- 


= { 0 } 


,..( 2 . 52 ) 

In the limit of y approaching zero, the cross receptances vanish 

i » e. , 

It 0 » i = 1,4 ...(2.53), 


Y ■+■ 0 
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Hence 


1 

1 

0 

0 


3y[(0) 

x,l 

e 

Np 

e 

0 

0 

< 

jyj (0) 

r 

0 

0 

1 

1 

1 

j 

y,3 

0 

0 

Np 

e 

■ Np 
e i 


yr4 


' ...(2.54) 

where the syston of equations has been decoupled into two 
sets of simultaneous equations which are identical to 
equations (2.44 and 2.47) obtained directly. As discussed 
in the following s ec tion , in ihe frequency regions where--- 
the coefficient matrix of equation (2.25) reduces to that of 
equation (2. 54)» one encoimters computatiox.al probl^s. 
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2.8 Results and Discussions 

2.8.1 Computations P erf ormed 

Using tlie analyses presented in the previous 
sections the following computations are perform eds 

(1) Curves for propagation constants versus frequency 

are obtained for = 0.6-n;, 6 = 0,3, y = 0 and y = 0.8. 

(2) The regions of principal primary instability are 
obtained with u^ = 0.6n; for the following values of the 
parameters in various cases. 

Case (a) ; By direct method for one span pipe with 
y = 0,8 and y = 0 in the range of 6 = O.i to 0,5. 

Case (b) t By the propagation constants approach 
using F(fi) for a one span pipe with the same values 
of the parameters as in (a). 

Case (c) ; By the F(n) approach for a two span pipe 
with the same values of the parameters as in (a). 
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table 2.1 

Tabulation of convergence rate with truncation 
order K at ditierent frequencies (u = O.671, 

Y = 0.8, 6 = 0.3h 


C onv er gene e Rat e ? ar am et er -q * j /i K= 4 


h 

20 

50 

35 

40 

45 

50 

14 

0 . 68 

2.14 

3.75 

7.18 

13.90 

26.30 

16 

0.074 

0. 22 

0.45 

1.041 

2.53 

5.32 

18 

0.007 

0.017 

0.029 

0.074 

0,19 

0.48 


Convergence Rate Parameter n* 

AE=4 


K 

50 


70 

90 

ICO 

110 

20 

0.04 


0.5 

9.0 

8.36 

- 

24 

0.00 


0.006 

0.117 

0.12 

2.14 

28 

0.00 


0.00 

0.00 

0.001 

0.04 

36 

n = 100 [ 

L 

p y 

(E + AK.) - 

s 

(K) I 

s 

-- 1 



L 

p y r. 

^ r , s 
s 

(E +AK) 1 

J 


where ^ is summation over p=x,y“, q=x,y; ri=0,l; s=0, 1; 

and 6^’^ (K' ) is the receptance value evaluated with truncation 
r , s 

order K = K’ . 
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2.8.2 Discussions 


The accuracy of the values for the receptances 
depends directly on the precision with which the truncated 
series solutions (equations (2.29a,h)) represent X(x) and 
y(x). The rate of convergence in the receptance values 
with increasing number of terms in the power series solution 
is presented in Table (2.1). At hi^er frequencies we 
require larger number of terms. The convergence rate is 
relatively insensitive to other system parameters. Since 
most of the computation is involved in the frequency range 
of s:2 = 0 to = 35 for all calculations the series is 
truncated at K = 16. For computing the propagation constants 
in the range q = 10 to n = 110, K = 32 is used. 


Figure (2.5) is a plot of the propagation constants 
and versus the non-dimensional frequency £2 for u^ = 0.6% 
Y = 0, 6 = 0.3. For the uncoupled case (y = 0) » propagation 

constants i = 1,4 are observed to occur in pairs such that 
= -^2? 12-2 "1^3 ^1 ^ ^ 3 * .noted that either 

purely real or purely imaginary as in the case of periodic 
beam [9], Considering the curve for \x^ (Fig. (2. 5)) it is seen 
that in the range £2 = 0 to £2=16. p^ is purely real and 


this frequency interval is identified as the attenuation 
band [9]. At £2 = 16.88 the real part goes to zero and over the 
interval 16.88< £2 < 43-3, P^ is purely imaginary and varies 
from % to 0. This zone is teimed as the 

the waves associated with lij propagate witagrtWjteni^^ 



^ 4 * -^*1 
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The propagation band is followed by another attenuation band 
and so on, alternately, as fi is increased. The propagation 
and attenuation bands for follow the same pattern as 
The plot for is displaced on the freq,uency axis relative 
to that of owing to the fact that the governing equations 

for X(x) and Y(x) even in the uncoupled case, are different 
(see equations (2.36 a,b)). 

Based on the discussion in section 2.7 the instability 
bounding frequencies may be obtained directly from the plot 
of p versus fi . Thus for N = 2, one set of bounding frequen- 
cies are obtained at Re (p^^) = 0 and Im (p^) = it, i = 1 and 
i = 3 i.e. the start of the respective propagation bands for 
p-j^ and p^. Similarly another set of frequencies for the second 
zone of instability are obtained with Re (p^) = 0, Im (p^) = 
i = 1 and i = 3 * In general for an N span pipe there will be 
N principal primary zones of instability clustered together 
within each set of propagation bands. It will be seen later, 
that the first zone of instability always corresponds closely 
to the frequencies marking the start of the first propagation 
band of p^ and p^ for all values of y and R. 

When y 0 equations ( 2.08a, b) are coupled and the 
curves for propagation constants also change. Fig. (2. 6) is a 
plot of P 2 _ and p^ versus « for the case of hi^ fluid mass 
ratio Y = 0.8 and u^ = 0.6n;, 6 = 0.3. ¥e observe that between 
the first and second propagation bands, and p^ are complex 
valued having both non-zero real and imaginary parts. Moreover, 



W ! 



Fig 2 7 Cutvflfs «^howinq u versus frequency In the first propasafron 
band 03 
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TABLE 2.2 

Values of frequencies Bounding Princioal Primary 
Instability Region for a Sin^e Span Pipe, Obtained 
by Direct Evaluation of frequency Determinant. 

Truncation order ; K = 16 

flow 7 elocity ; Uq= 0.6% 

fluid Mass Ratio ; y = 0 (uncoupled case) 


Y = 0.8 (coupled case) 
fRE QUMO IBS 


Excitation | 
Parameter 

6 

i — - - * 

Y = 

i 

0 

Y = 

0. 8 

Lower 

Upper 

lower 

Upper 

0.01 1 

I ! 

1 - , . 

15.747 

15.836 

15.347 

15.506 

I 

i 0*10 

15.318 

16.208 

14. 614 

16. 203 

8. 20 

14.782 

16.572 

! 

13.766 

t - - 

16.945 

j0.30 

1, 

14.177 

16.887 ' 

12.881 

17.650 

1 

'0.40 

i 

1 

13.492 

17.155 

11.955 

18. 319 

^ 0.50 

f 

f 

; 12.716 

{ 

\ 

17.378 

10.982 

18. 950 

1 

i 
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in this frequency range (i.e. t-i^ and occur as 

conjugate pairs). As in the earlier case the set of |j,^, 
i = 1,4 remains closed under the unary operation (i.e. 

M-l = -1^2 and at all frequencies. 


Figure (2.7) is a plot of over the first 

propagation hand for y = 0.8 and y = 0 with u^ = 0.6%, 6 = 0.3. 
On comparing curves of p^ and p^ for y = 0 and y = 0. 8 we 
observe that the effect of increased coupling is to magnify 
the relative shift between the curves of and p^. Increased 
separation between p^ and p^ plots indicates broadening of the 
instability zones as will be evident shortly. 


The validity of the graphical approach in deteimiining 
the instability bounds is restricted only to the uncoupled 
case (i.e. y = 0). To obtain the instability bounds for the 
coupled case (y 0) we take resort to equation (2.25). 


The frequency bounds for a single span pipe are first 
obtained for both the coupled and uncoupled case by the direct 
method mentioned in section 2.3* Table (2.2) is a list of these 
results for u^ = 0.67c, y = 0.8 and y = 0, N = 1. These frequen- 
cies closely correspond to the start of the propagation band 
of p^ and \iy For the particular case of 6 = 0.3 this may be 
compared by tallying the pair of frequencies in Table (2.2) 
with Fig. (2.7). The results obtained by this method are 
in excellent agreement with those obtained in reference [7] 
where the Hitz-Galerkin approach was used. 
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MLE 2. 3 

Values of Frequencies Bounding Principal Primary 
Instability Region for a Single Span Pipe, Obtained 
by Using Propagation Constants and the Determinant F(Q.). 


Trane at ion order s K = 16 

Flow Velocity ; 0,&k 

Fluid Mass Ratio : y = 0 (uncoupled case) 




Y = 

FREQUEIIC lES 

0.8 (coupled 

case) 

! 

Excitation 

Parameter 

6 

■y 

= 0 

Y = 0 

.8 j 

Lower 

Upper 

Lower 

Upper j 

■■ 0.01 

15.75 

18.84 

15.35 

15.51 

0.10 

15.32 

16.21 

14.61 

16.23 

0.23 

14.78 

15.57 

13.77 

16.95 

0.30 



1 2. 881 

17.650 

1 

; 0.40 

1 

13.49 

17.16 

12.00 

18.32 

1 

0.50 

12.72 

17.38 

10.98 

18.95 
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In tile frequency deteim inant approach, the bounding 
frequencies are obtained by searching for the zero crossings 
of P (n) (equation (2.25)). lo this end P(n) is plotted against 
Jli. and a sample plot is shown in Pig, (2.8) for u^ = 0,6%, 

6 = 0, 5, H = 1, Y = 0. 8, It is observed that at any frequency^ 
Xl the function P(n) is either purely real orpur«ly imaginary. 
The function exhibits zero crossings of the following types: 

(1) P(n) remains purely real or purely imaginary in the 
vicinity of n . where 

J 

P( .) = 0. 

(2) P(n) changes phase from real to imaginary or vice-versa 
as it changes sign across the point = n . where 

J 

P(fi-) = 0. 

At the frequencies, fi = 12,881 and 17,650 the function 
exhibits zero crossings of the first kind and at n = 15-02 and 
18,109 zero crossings are of the second kind. 

On comparing these frequencies x-irith those obtained by 
the direct method for the same set of parameters (Table (2.2)), 
we infer that only the zero crossings of the first kind corres- 
pond to the bounding frequencies. The zeroes of the second 
kind exist due to some numerical problem discussed later in 
this section. Table (2.3) lists the results obtained by this 
method for N=l, u^ = 0,6 tc, y=0f 0,8. These results tally 
exactly with the values in Table (2.2). 
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^’igure (2.9) is a representative sketch of ?(«) versus 
2 with Uq = 0, 6jr, y = 0.8, 6 = 0,3. The first pair 
of aero crossings (first kind) in this figure at Si = 12.90, and 
^ = 17.53 yield the bounding frequencies for the first zone of 
instability. With reference to Fig. (2.7), these frequencies 
closely correspond to the start of the propagation bands for 
and (y = 0,8), Thus for the first zone even for y = 0,8, 
the results are close to those obtained by the graphical approach 

The second set of zero crossings in Fig. (2,9) at 
n = 25.25 and si = 29.71 correspond to the second zone of instabl 
lity. At these frequencies, 

Re (p^) = Ee (p^) = 0 

but 

(P-l) ¥ %/2 ¥ Im (p^) 

which indicates that the graphical approach cannot be used to 
evaluate the second (or higjier) zone of instability for y ^ 0. 

From this discussion we may conclude that the first 
zone of instability bounds always occur close to the start of 
the respective propagation bands for p^ and p^. This is 
regardless of the number of spans in the structure or the degree 
of coupling between X(x) and Y(x), as concluded in reference [12] 
F or y 0 case, the hi^er zones however cannot be obtained by 
the graphical approach as concluded in reference [12], These | 
zones may be evaluated accurately only by using equation (2.25), 




It is shovna in section 2.7 that the graphical approach 
applies when the matrix in equation (2.52) becomes decoupled. 

The ft ct that at the bounds of the first instability zone, the 
graphical approach and the frequency determinant approach, giv e 
identical results suggests, that at these frequencies equation 
(2.52; is effectively decoupled. This implies that even for 
Y = 0.8, r^ = 0, (i = 1^4) near the start of the propagation 
bands of V-y In other words cross receptances are very 

small in comparison to the direct receptances at these two 
iroquencies. However the importance of the cross receptances 
while calculating the propagation constants is easily seen from 
the considerable change in the pair of curves (p^,* 
two cases y = 0.8 and y = 0 as shown in Pig. (2*7)« 

In the vicinity of the boundaries of the first zone, 
since r!^ 's approach zero, irrespective of the value of y, the 
determinant F ( q) becomes ill conditioned with one or more rows 
having very small numbers in comparison to the other rows. This 
results in the spurious zeroes (zero crossings of the second 
kind) motioned earlier. No spurious zeroes were encountered 
near the bounds for the higher instability zones. Also these 
results are significantly different from those obtained by 
graphical approach, indicating the increased importance of the 
cross receptance for these zones. 

Figure (2.10)is a plot of the instability zones in the 
n - 6, plane for 1 = 1,2, y - 0^ 0.8 and u^ = 0,6it. 
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TAELS 2.4 

Values^of frequencies boimding the principal primary 
instability regions associated with the first two modes 
of a tvro span pipe (using propagation constants and 
determinant E(a)), K = 16, u^ = 0.6tc 

Frequencies (Uncoupled case y = 0) 


Excitation 

1 I mode I 

II 

mode 

P arameter 

6 

L ower 

Upper 

L ower 

Upper 

0.01 

15.747 

15.836 

28.023 

28. 082 

0.10 

15. 318 

16. 208 

27.743 

28 . 330 

0.20 

14.782 

16.572 

27.399 

28.576 

0.30 

14.177 

16. 887 

27.020 

28.791 

0.40 

13.493 

17.455 

26 , 6 O 4 

28.976 

0.50 

12.716 

17.378 

25.149 

29.132 

Frequencies 

0 . (coupled case y = 

0.8) 

Excitation 

I mode 

II 

mode 

Parameter 

6 

lower 

Upper 

1 lower 

Upper 

0.01 

15.348 

15.598 

27 . 463 

27.612 

0.10 

14.625 

16.183 

26.773 

28.289 

0.20 

13.781 

16.887 

26.006 

29.01? 

0.30 

12. 897 

17.528 

25.254 

29.714 

0.40 

11. 969 

18. 083 

24-540 

30.374 

0.50 

10.994 



18.516 

23.905 

30.993 
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In the range of frequencies considered we ohserve 
one and two instability zones for N = 1 and 2 respectively. 

From these curves we conclude that the width of the instabi- 
lity zones for both y = Y = 0 increases with 6» On 

increasing y, from 0 to 0.8 the instability zones become 
broader and this difference is seen to be greater for higher 
6’s, In the low 6 region (0 - 0.1) it is observed that 
increasing y shifts the zone to lower frequencies. These 
conclusions are at a variance from those obtained in [5]and[l2} 
where it was observed that the effect of increasing y is to 
depress the instability zones uniformly to lower frequencies 
and that this change is less for hi^er values of 6. ®ie 
errors in reference [12] arise due to 

i) Improper analysis at the stage of obtaining the 
propagation constants. 

ii) Improper equations in deriving the end momoit 
conditions, 

iii) Use of a different governing equation based on 

reference [6] which itself is erroneous as shown 
in refer enpe [7]. 

Numerical values for the two span pipe for = 0*67 e, 
y =5 0, 0,8 are presented in Table {2,4)* 
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CHAPTER - III 

CQtJTROL OF' parametric INSTABILITIES BY DYNAMIC ABSORBERS 
3.1 Introduction 

Dynamic absorbers are used to control the excessive 
vibrations in mechanical systems. By suitably varying the 
absorber parameters, the system impedance can be changed 
accordingly. Unlike dampers which reduce the system response 
by dissipating energy, an undamped absorber diminishes the 
response by force cancellation, A dynamic absorber normally 
consists of a lumped mass that is attached to the main (or 
primary) systen by a damped resilient element. 

In the present chapter, a viscously damped dynamic 
absorber has been, tried for a different use from its conven- 
tionalrole of r evincing the resonant motion. The aim of the 
present study is to determine whether the dynamic absorber 
can effectively control the regions of parametric instabilities 
of a pipe. 

As stated in Chapter - II in most of the cases, the 
principal region of instability is of primary importance. The 
fluid mass ratiO> parameter, y bas been retained since it has 
been found to ® the primary zones significantly 

(Section 2.8), maintain the periodicity in the overall 

structxjre, the n-'bsorber need be attached at identical location 
in each span, present work, the point of attachment of 
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the absorber is taken as the mid-pohat of each bay. Computa-* 
tions have been perfomed for U 2 only, thou^ the 
method used is equally valid for any number of spans. The 
elegance of the method lies in the fact that only the matrix 
used in Chapter-Il for receptance calculation need be changed. 
The rest of the computation is on identical lanes with the 
case without the absorber. 

3 » 2 Equations of Motion and Receptances 

Consider a pipe of M-spans, each having a length I 
rooting on line supports (as discussed in chapter-II). a 
viscously damped absorber of mass ma, spring constant k and 

3 . 

viscous damping coefficient is attached to the mid-point 
of each span (Pig. ( 5.1a)), The periodic element for this 
structure is shown in Pig. ( 5 . 1 b). let w, (x,t) and ^ 7 ^(x,t) 
refer to the transverse displacement of the left (x = 0 to I/2) 
iialf and the rigiit half (x = 1/2 to 1) respectively, of the pipe 
ol em ent . 

¥e shall, proceed first to derive the boundary 
conditions at x = 1/2 imposed by the addition of the absorber. 
(Pile displacement, slope and mcoent continuity relations 

I 

across the absorber are, 


...( 2 . 55 a) 
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m* 

9X 





...(2.5513) 


^2 - 
9 Wt 


3X 


( x) 


^ 

-2 '■ 2 * 
8X ^ 


t;) 


...(2.55c) 


Following u non-dim ensionalising scheme as in section 2.2 
we may rewrite equations (2.55a,b and c) as 


Wj. ( |, t) = Wjj ( i, t) 


"l ^ 2’ 


( 2» "t) 


...(2. 56a) 
...(2.56h) 


“i' ( 2* 


( 2> 


...( 2 , 56 c) 


where and are the non-dimensional equivalents of ^ 
and w^, respectively. 

For the K = 1 approximation (section 2.3) we may write 


Wt 


sin ( i S2 t ) + cos ( -^ a t) 


...(2.57a) 


w. 


R " "R 


Xj , Sin ( i n t ) + Tjj oos ( i St) ...(2.57b) 


Using equations 


(2.56a,h,c) and (2. 57a,b) we obtain 


Xi ( I ) = Xi ( I ) 


..,(2. 58a) 


. . . (2.58b) 
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H' 

(1) 

= ^r' ( 1 ) 

. . . (2.58c) 


( |) 

^ ^ i ' 

. ..(2.58d) 

Y I 

1 

(|) 

= Y' { ~ ^ 

...(2.58e) 

Y ' • 

^1 

( 1 ) 

= Y» » f i. \ 

R ^ 2 ^ 

... (2.58f) 


The shear force balance at 3C = “ yields the fourth coupling 
relation between Wj^ and w^^. let Vp and be the shear forces 
in the pipe at the ri£^t and left of the absorber location 
(at i = 1/2). Then 

•— 3 W-r- T 

= El -^ 3 ^ ...(2.59a) 

3X 

_ 3 ^ Wp T 

¥p = t;) ...(2.59b) 

Let y be the independent coordinate describing the transverse 
displacement of the absorber mass (Fig. (3.1b)). 

Balancing the shear forces at 5 = ^, together with the 
equation of motion for m^,we obtain 

. . . ( 2. 60a) 

and ^£“"^ 1 , = “ “a ...(2.60b) 
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Jb'-L Oil i iici 'G ^ 


' 7 v 5 <Dbtair] 


I and y tran eciuatlons {2,5Sa,-D) and (2.6Ca,’o) 


iiX 


cr 






Ox 


3 *■ °‘a + C_ iL +■ t, “ ■) ^-r. 3^ 


a '"r + J-a J = BI [in 

3 s 




3 -o Wt 4- K “r_ 1 

-r-^ d. 1 j a 1 ‘a I-' 


'"a \ \ ] 


ax 


1/2 


( 2. Gl) 


• . * .r e ( ,< ocnoteii derivative 
i.'<-. r* t-‘kj'->.i,nA, aijovo equation in 


with respect to time, t. 


a non-diuieiiGional foim,- we obtairs 


t '^K + 'a H. , ^.4 5 

3x a a 


W-r 

jj 


+ '='a r [ 


1 + \ at 2 = 1/2 




'fiio 'e ( ) - 3 “^- , non~diaeneional time derix^ativ 

C 


a \K m~~ 
a a a 


, the damp in £ factor. 


O’ 


a 


j,2 i<.,. (m^+m ) 1/2 

2 ^ m El f the xuixing ratio^ 


n; 


m 




a. ~ Jm., + in )l ' absorber mass ratio. 

*P 



es'pressioas for 


Substituting the 

b) in 


equation; 

i above equation we get, 


• 



r* 

t 





1 

o 

H 


i ..'"1 

1 X,- 1 

1 1., 

t 

o 

VO 

"2i 

t tt ! 

i""! : 


1 

i 

1 

>- 

1 

+ 1 

j ' 

' ‘ f 

, ? 

'''2 

n . 


i 1 

ui 


i Q 

1 

“"ij I 

! It. 

i Y ■ 

1 ^ 


D. 


D, 


D- 


'1 

X = 1/2, 








r ~ 


0 


• • ( 2» 6^ j 


wh or e 


C 


1 


[ 3. 


f2c 


2 4 ] » 3)^ 

4o-^ TC ^ -*■ 


D, 


fj 

3'a 4 


2 


”a 8 


Equations ( 2, 5Ga., b, c,d, e,f) and (2,63) are the eight equations 
v.'n Lcri couple the solutions in the right and left haif of the 
p ip e, 

Tb.e conditions at the ends jc = 0, x = 1 with reference 
to section 2.6 and i'ig. (3.1b) laay be directly written as. 



(0) = 0 

, (1) = 0 

. . . ( 2. 64a) 


(0) = 0 

, . (1) 0 

, . , (2.64b) 
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H' ( 0 ) = , XI' ( 1 ) 

Ji' ( 0 ) = - k | , !■. ( 1 ) 




. . . ( 2. 65a) 

% 


. . . ( 2, 65t ) 


The governing equations that apply to the left and. ri^t 

half of the pipe element are respectively (see equation (2.08a,b)). 



+ u: 


(1 - 6 + 




1 

4 


2; 

1 


YUq6 


[(l~x) 




0 . , . ( 2. 66a) 




dx 




d^Y, 


dx‘ 


7 I 
4 






[(l-:5c) 


d^X, 


dx' 


o dZ 

- + 1> di^ ] = 


0 ..,(2. 66b) 


where the subscript q = L with 0 <_ x £ 1/2 

= R with 1/2 < X £ 1. 


Equation (2.66a) and (2.66b) are four, fourth order coupled 
differential equations and together with the sixteen constraints 
given by equations ( 2. 58a,b,c,d, e, f ; 2«63j 2. 64a,bj 2.65a,b) 
form a determinate system of equations for the four imknowns 

f Yj^ , 2^ f . 
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The above system of equations may he solved by express- 
ing each unknown as a truncated power series. The receptances 
are obtained in the same manner as outlined in section 2.6. 

3-3 Results and Discussions 

3.3*1 Computations p erf ormed 

Dsing the analyses presented in the previous sections, 
mjimerical results, listed b^ow have been obtained. As in 
chapter II, the instability regions are presented as plots of 
non-dimensional frequency Q. versus the excitation parameter 6. 
All results pertain to the following pipe parameters : 

= 0.6 Ti, R = 2, Y = 

2 

(1) Determination of instability regions with =0,2, = 1 

and damping factors ? = 0,01, 0.05, 0.1. 

(2) Instability regions are obtained to study the effect of 

changing the tuning ratio (o^ = 0,8, 1, 1.2) of a 
li^tly damped absorber ( = 0.01, = 0.2). 

Results presented in chapter II are used as a check, for the case 

r = 0, in the absorber formulation. 

a ’ 

5,3.2 Discussions 

Convergence in the recc^tance values in frequency 
range of interest was obtained by truncating each powers series 
of \ after 16 terms. 



irj 32 Rigions principal primary instability 

* — without absorber Uo‘0'6n 2 0*6 
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F ig. 3*3 Effect of damping factor on the first pair of principal 
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Fig. 3’4 Effect of tun ir^ ratio on the first pair of princip-f 
primary instability zones (absorber case) 

with absorber uo=06inf,N=2 ,T=O-0 

f;=o-2-Sa=o-oi .ol-o-s.vo.i; 

— Without absorber 00^-06 n jN-2 , T'-O F 



Fig 3-5 Effect of tuning ratio on second pair of principal 
primary instability zones (absorber ease) 

— — with absorber Uo=06Tr,N-2,T=0*8 

r;=0-2,'Sa=0'01 ,0^=06 ,10, V2 
Without absorber iu-0'6ti»N=2, y=0-8 
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Fig. (3»2) is a plot of the instability zones for 
N = 2 , U- = 0.6 Tt, Y = 0.8 and r = 0 , 2 , = 1.0 and ? = 0 . 01 . 

V 3 Q, S' 

In the same figure for comparison the instability zones for the 
same primary system parameters, without the absorber are also 
shown. 

¥e observe that each instability zone (in the non- 
absorber case) splits into two, so that for a two span pipe we 
havG four instability zones in the region of investigation, 
ihese zones are referred to as A 2 > 82 * for convenience, 

as shown in Fig. (3.3). Ihese zones are narrower in width in 
comparison to the instability zone without the absorber (marked 
A and E in the same figure). It is noted that A 2 anti B 2 s^e 
narrower compared to A^ and The zones A^, B^^ are further 

apart on the frequency axis in comparison to A and B. 

Besides the repositioning of the instability zones 
along the frequency axis on addition of a li^tly damped absorber 
it is observed that the zones A 2 » ^ 2 ’ ®1 shifted along tiie 
6 axis too. This is ascribed to the eif ect oi damping (present 
in the absorber ) [l£j . Thus with each instability zone there is 
on assosiated critical excitation parameter 6 * (e.g. with 

A ) such that the loss of dynamic stability (on account of a 
particular zone) occurs only when 6 > 6 . 

At C = 0.01 (which is a case of low damping) we 
a 

observe that 6 * = 0. In the range of parameters investigated 

•^1 

it appears that gonerally 
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and 6* ^ 6* 

It is seen that, hi^er modes are significantly shifted along 
the 6 axis. To observe the effect of damping (particularly for 
zone is increased to 0.05 and 0.1. It is seen that 

on changing frcan 0,1 to 0.05, the instability zones A 2 » 

B 2 are swept off to the ri^t along the 6 axis so that 

(6f , 6* 6* ) > 0.8. Figure (5.3) is a plot of the I mcde 

•‘=2 ^1 . 
instability zones for various values of From this figure 

we observe that zone A^^ can also be successfully controlled 

such that 6* > 0. The instability zone A^ is a relatively 

’hard zone' and is affected to a lesser extent by the absorber 

damping. Thus the qualitative effect of increased dampitig is 

to shift the zones towards the ri^t and this effect is more 

pronounced for hi^er modes of instability. At this stage it 

is remarked, that the critical excitation parameter 6 , is also 

a fimction of the other absorber parameters and in general, 

changes on varying and as well. Ihe effect of toning 

ratio, on the instability zones is shorn in Tig. (3-4) and 

( 35 ). These plots correspcnd to the first and second pair 

of modes respectively with ^ q . "" 

1,0 and 1.2. For these values of parameters the effect of 
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2 

changing on zcne is n eligible compared to that on ^ 2 * 

p 

On increasing o^, A 2 is observed to shift to hi^er frequency 
and that 6. also increases. The effect of changing cr„ is 

Aq O' 

c 2 

Opposite on and B 2 (yig. 3- 5) ) . On increasing from 0,8 

to 1.2, zones Bp and B 2 shift to lover frequency and 
-K 3E 

6. and 5? decrease in magnitude. 

"1 ■“2 

This suggests that with c fixed there exists an optimum 
value so that both 6^ and 6^ , i = 1,2 are a maximum, 

The effect of change in the absorber parameters is thus 

albeit, 

seen to significantly alter that instability zones^in a complex 
manner. It may be concluded from these investigations that 
the effect of damping is more pronounced for instabilities 
in the higher modes. The optimisation of tuning ratio may 
be necessary in the range of parameters where affects 
and B 2 _ significantly. 

In actual conditions, the fluctuations in the flow are 
not harmonic. They may be described as a band limited 
stochastic process. Once the power spectral density of the ; 
fluctuations is obtained we can define the region of inter sst 
in n - 6 space that is to be made entirely stable. By prescribin 
a suitable objective function, for e. g. , (i) maximize 6 , 

(ii) minimize width of instability regions, (iii) maximize 
frequency separation of instability zones, etc,, we can I 

suitably optimize the absorber parameters for best results. 
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CHAPTER - IV 

CQMCLUSIONS 


4.1 Concl-usions 

Prom the results presented in Chapter 2 and 3 the 
following conclusions are drawn. 

The parametric instabilities of a periodically- 
supported pipe conveying a piiLsatile flow can be convenien- 
tly analysed by the propagation constants approach. This is 
mentioned particularly, in view of the appearance of the space 
dependent term and the flow direction dependent terms -vdiich 
destroy the directional equivalffipce in the governing equation- 

The coupling effect of the fluid mass ratio parameter 
has been incorporated in the analysis in a systematic fashion 
(by introducing cross rec^tances) and the hitherto unknown 
effects have been discussed. The propagation constants have 
been shown to be sensitive to the fluid mass ratio and it is 
shown that for high values of this parameter the coupling 
effect must be retained at all stages of analysis. For low 
degree of coupling it is shown that this approach is identical 
to the graphical approach available in the literature. The 
effect of increasing the fluid mass ratio is seen to broaden 
the instability zones and that this effect is more pronounced 
for higgler values of excitation parameters. 
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It is Seen, that irrespective of the degree of 
coupling, the frequency pair bounding the first zone of 
instability corresponds to the start of tbs propagation bands 
for th e respective propagation constants. The frequencies 
bounding the higher zones of instability do not bear a sinple 
correlation with the plot of the propagation constants versus 
the non-dimensional frequency n (especially for the case of 
high fluid mass ratio). F or a low degree of coupling, however, 
the graphical technique may usefully be applied. 

The instability zones are markedly affected by the 
attachment of dynamic absorbers at the mid-point of each span. 
The methodology in the analysis is unchanged on addition of 
absorbers to the system, if the periodicity is retained. In 
the same range of frequency, the nimaber of instability regions 
with absorbers are more, than those without it. However, the 
width of these regions are much smaller than the original 
regions and they can be controlled by varying the absorber 
parameters. With a fixed absorber mass ratio parameter, the 
first zone of instability is affected to a lesser extent by 
the change of the timing ratio and the damping factor. The 
critical excitation parameter, in general, increases with 
hi^er values of the damping factor. The trend of the effects 
of the tuning ratio, on the different zones of instability, 
suggests the existence of an optimum value. Hence depending 
upon specific requirements, the optimum values of these 
narameters can be determined. Due to the breakdown of 



directional parity in the system (on accomt of flow) the 
best location for the absorber needs to be ascertained in a 
systematic way. 

The method developed in this work can readily be 
extended for determining the secondary and hi^er order 
(K = instability regions as well. Moreover, pnlike 

the Ritz-G-alerkin procedure [7] the primary instability 
regions associated with all the modes are obtained simulta- 
neously. 

4.2 Recommendations for Future Work 

In view of the simplifications and the results 
obtained in the present thesis, the following directions 
are suggested for further research. 

i) Incorporation of the effect of gravity, pressuri- 
sation, axial tension and contraction which have 
not been explored in the present work can be 
undertaken. 

ii) Parametric regions (with pulsatile flow) that 

trigger combinational resonances can be investigated, 

iii) Optimisation techniques for absorber location and 

absorber parameters for control of instability zones 
can be investigated. 

iv) The problem of parametric instabilities of pipes with 
random fluctuations in the fluid velocity can be 
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